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Abstract 

The KelvinHelmholtz (KH) instability is a fundamental wave instability that is frequently observed in all kinds 
of shear layer (jets, wakes, atmospheric air currents etc). The study of KH-instability, coherent flow structures has 
a major impact in understanding the fundamentals of fluid dynamics. Therefore there is a need for methods that 
can identify and analyse these structures. In this Final assignment, we use machine-learning methods such as Proper 
Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) to analyse the coherent flow structures. 
We used a 2D co-axial jet as our data, with Reynolds number corresponding to Re: 10,000. Results for POD modes 
and DMD modes are discussed and compared. 


I. Introduction 


The KelvinHelmholtz (KH) instability can occur when there is passive scalar shear in a single continuous 
fluid, or where there is a passive scalar difference across the interface between two fluids |[TJ. To gain 
deeper understanding on KH coherent flow strcutures, machine-learning methods can be utilised. This field 
allows computers to adopt behaviors based on training data. These methods recognize complex patterns 
and makes intelligent decisions based on data Q. These techniques include: reduced order models or 
dimensionality reduction methods Q, statistical, and machine vision methods. Dimensionality reduction is 
a technique used to find a reduced order model on a given data. Such a technique includes feature selection 
and feature extraction methods. Feature selection is based on selecting a subset of variables which best 
define the data, whereas feature extraction transforms the data from high-dimensional space to a space 
of fewer dimensions Q. The data transformation may be linear, as in Proper Orthogonal Decomposition 
(POD) Q 0 ||6j and may be non-linear such as in Dynamic Mode Decomposition (DMD) Q @ |[6|. 


From the literature, it is expected that POD and DMD would be good methods for analyzing the flow 
structures. For example, Perrin et. al Q used POD to obtain phase averaged turbulence properties for flow 
past a cylinder. In highly turbulent flows, the coherent flow structures are difficult to identify due to the 
combination of organized and chaotic fluctuating motions. Using POD analysis it is shown in Q, that 
von Karman vortices can be reproduced within the first few modes. POD has been used as a tool for the 
comparison of Particle Image Velocimetr 5 |^PIV) and Light Eddy Simulation(LES) data in i fT^ and it is 
also shown that POD modes have a good qualitative agreement between PIV and LES. A paper by Schmidt 
et. al (H). used DMD to a sequence of flow images of a slow jet entering quiescent fluid showcased the 
detection of dynamically relevant coherent structures that play an important role in characterizing the fluid 
behaviour over processed time interval. To identify the KH coherent flow structures, we study 2 cases of 
co-axial jet: 


• When both the jets are close to each other (CASE = ^). 


'http://en.wikipedia.org/wiki/Paiticle_iniage_velocimetry 




• When both the jets are independent, far a part from each other (CASE = f). 

where L is the domain side length i.e 2% or 256 pixels, as the grid resolution N is 2® pixels. Both the 
cases of co-axial jet are studied at Reynolds number corresponding to Re: 10,000. The inverse momentum 
thickness or the steepness parameter B = 10.5, the jet diameter Vo is kept constant at The U,V and PS 
for a co-axial jet are defined as follows: 

Umax = 0.1; Vmax=Umax/30; ro = L/20; 
yl= L/2 - CASE; y2= L/2 + CASE; 

U = Umax*0.5*(1-tanh(B*(abs(Y-yl)./ro - ro./abs(Y-yl)))); 

U = U-i-Umax*0.5* (1-tanh (B* (abs (Y-y2) ./ro - ro ./abs (Y-y2) ) ) ) ; 

U = Umax*U/max(max(U)) ; 

V = Vmax*(rand(N,N)-0.5) ; 

PS = U/max(max(U)) ; 

Figure shows the mean profile for the two cases investigated in this report. The interactions between 
the jets in case 1 starts at step 40, whereas in case 2 it is at time step 100. 
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Fig. 1. Mean Profiles of case 1 (left) and case 2 ( right) passive scalar fields 


The objectives of the paper is as follows: 

• POD and DMD is implemented with Matlab and used to analyse LES of subsonic jet. 

• The potential of POD and DMD in the ICE applications is pointed out. 

• Results for POD and DMD are discussed and compared. 

The organisation of this report is as follows: In section 2 we study the computational methods for POD 
and DMD. Section 3 deals with the experiments and implemention details. Finally, conclusions are drawn 
and further research directions are discussed in section 4. 

II. Computational Methods 
In this section, the POD and DMD methods are discussed. 





A. POD (Method of Snapshots) 

POD was first introduced in the field of Computational Fluid Dynamics (CFD) by Lumley [ [T2| |. POD is 
known by various names like KLD, Principle Componant Analysis (PCA) and Singular Value Decompostion 
(SVD). This method aims at representing a 2D or 3D flow field in terms of orthogonal modes which 
correspond to the energitic, coherent flow motions. Hence, POD may be used to extract certain derterministic 
features from turbulent data. POD theory is best described in Q. The present day analysis uses method 
of snapshots introduced by Sirovich (T^. Here LES data at a particular interval of time, interpolated on a 
uniform grid is considered to be one snapshot. This method was introduced to reduce the POD computations. 
To compute the POD using the original method requires solving nxn eigenvalue problem. The main problem 
is the calculation of auto-covariance matrix R. Method of snapshots proposes that auto-covariance matrix 
can be approximated by a summation of ’snapshots’. Here x is snapshot matrix. 
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The snapshots are assumed to be distanced of spatial distance greater than the correlation distance. Nowadays 
this method has become extremely popular [ fT4| | Q [ fT5] |. Its use in certain flows are questionable however, 
due to its assumption of the snapshots being uncorrelated. This discussion is briefly explained in 0 


B. Dynamic Mode Decomposition (DMD) 

The mathematics underlying the extraction of dynamic information from time-resolved snapshots of LES 
data is closely related to the idea underlying the Arnoldi algorithm. Let Vj denote each instantaneous flow 
field. A sequence of N is writen as : 

Vf = [t;ir;2V3t'4 ■■■ vn ] (2) 

A linear mapping from one snapshot to another is assummed. 

Vf = [viAviA\iA\ ■ ■ ■ A^-^i] (3) 

This can further more assumed to be a constant, so that the mapping A can be represented as: 

Vj+i = Avj (4) 

By the linear combination of available data fields, we have a standard Arlondi iteration problem 0 

Av^-'^ ss Jif-'S (5) 

where S is a companion matrix that simply shifts the snapshots 1 through N — 1 and approximates the 
last snapshot V by a linear combination of previous N — 1 snapshots. Hence this procedure will result in 
the low dimensional system matrix S. We solve the S matrix problem using eigenvalue analysis and obtain 
the eigenvalues. It is known that eigenvalues of S, approximate some of the eigenvalues of the full system 
A. The associated eigenvectors of S provide the coefficients of the linear combination that is necessary to 
express the modal structure in the snapshot basis. S matrix is calculated as follows: 

S = R-^Q*vn (6) 

where Q* is the complex conjugate transpose of Q from the QR-decomposition of 


III. Experiments and Implementation 

Here we discuss the Experimental and implementation details of POD and DMD on 2D co-axial jet data. 



Fig. 2. case 1 (left) and case 2 (right). 


A. Pre-Processing 

Figure shows 2D co-axial jet for both the cases. A total of 30 snapshots were considered with a time 
interval of 5 steps. 

Consider an example of a smoke in the wind flow, the concentration of the smoke C at a particluar time 
t, is given by (7 = C{P, t). As part of experiments, snapshot matrix X is constructed at first. Let P be the 
scalar concentration component: 

C = [P^p2p3 ...pN^ = 

B. Implementation of POD 

The fluctuating passive scalar matrix U is calculated by substracting the mean from individual snapshots. 
Then the autocovariance matrix is computed as: 

C = U^U (8) 



The eigenvalue problem for the matrix reads as follows: 

CA^ = AM* (9) 

The eigenvectors are arranged according to the decreasing order of eigenvalues reflecting the energies in 
different POD modes. 


> • • • = 0 . 


Using the ordered eigenvectors the POD modes are constructed. 
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C. Implementation of DMD 

Calculation of fluctuating passive scalar matrix is not needed. But the snapshot matrix is divided into two 
parts. 


Vf ^ = [viV2Vz • • • Vn-l] 

(12) 

Vf = [V 2 V 3 V 4 ■■■ Vn] 

(13) 

QR decomposition in economy mode is performed as : 


[Q,R] = qriVr\0). 

(14) 

Companion matrix S is calculated as : 

S = R-^Q*Vf. 

(15) 

The eigenvalue analysis is computed on S matrix 


[X,D] = ezg{S) 

(16) 

The dynamic mode spectrum is computed as : 


Xj = log(Djj)/St. 

(17) 

St, is the time interval between the snapshots. The Dynamic modes can be 

computed as follows: 

DM, = 

(18) 

where X is the original snapshot matrix. Dynamic mode decomposition 

contains not only information 


about coherent structures, but also about their temporal evolution. Since at no stage of the algorithm the 
system matrix A is needed, various extensions and attractive features of the algorithm should be noted. No 
speciflc spatial arrangement of the sampled data is assumed, and the processing of subdomain data, i.e., data 
extracted only in a small region of the complete flow domain, as well as the processing of unstructured data 
is possible. This feature allows to probe a subregion of the flow as to the presence of a specific mechanism. 

IV. Results and Discussion 

The POD modes are the optimal decomposition for the flow and capture large scale structures thereby 
providing information on the large scale behaviour. Total kinetic energy is contained within the first few 
POD modes. Generally 95% of the total intensity fluctuations is used and the average flow field is described 
in the first n POD modes. Dynamic modes represent the perturbation dynamics and capture the characteristic 
pattern located near the shear layer. DMD spectrum quantatively describes the jet behaviour. Higher the 
dynamic modes higher the is the frequency. 

From the Figure it is observed from POD modes, that there is a formation of symmetric structures of 
inverse Kelvin Helmholz instability due to interactions of both the jets. But from Dynamic modes, the 
inverse Kelvin Helmholz instability is clearly visible. The second dynamic mode shows a characteristic 
pattern located near the shear layer, which represents the roll-up of the symmetric vortex sheet in Figure 
The first and third mode depict small scale structures near the nozzle exit region. 

From the Figure it is observed that there is the formation of the Kelvin Helmoltz instability in both the 
jets. From Figure the time coefficients of POD modes 3,4 show a phase difference of 5 time steps. In 
the other words, the POD modes 3,4 are strongly correlated to each other. This phenomena would further 
give rise to study helical and columnar features of instabilities. 

POD and DMD methods can be utilised for analysing the LES data. The power of PODs fast convergence 
allows for the large scale structures to be isolated from the small scale structures in the turbulence. This 






Fig. 3. Comparison of first 5 DMD modes (top) and POD modes (bottom) for case 1 passive scalar field. 



Fig. 4. Comparison of first 5 DMD modes (top) and POD modes (bottom) for case 2 passive scalar field. 



Fig. 5. Time coefficients of case l(left) and case 2( right) passive scalar fields 


would help in analysing the flow field in different ways. From the Figure and Figure POD and DMD 
methods are shown as a powerful numerical tools, for use in jet flows. The ability to maximize the intensity 
fluctuations of the flow with a miminal number of modes shows PODs strength in the analysis of the 
coherent flow structures and reduced order modelling. The DMD method has clear advantages over POD, 
as it strives for a representation of the dominant flow features with in a temporal orthogonal framework, 
while POD is based on a spatial orthogonal framework. 

From the Figure the eigenvalues of S represent the mapping between subsequent snapshots: unstable 
eigenvalues are given by a modulus greater than one (i.e. |A| > 1, are located outside the unit disk); stable 
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Fig. 6. Eigenvalues of DMD for case 1 passive scalar field. 


eigenvalues have a modulus equal to one (i.e. |A| = 1, can be found on the unit disk); unstable eigenvalues 
have a modulus less than one (i.e. |A| < 1, can be found inside the unit disk);. For applications in fluid 
dynamics, it is common to transform the eigenvalues of S using a logarithmic mapping, after which the 
unstable (stable) eigenvalues have a positive (negative) real part. The procedural steps for computing the 
dynamic mode decomposition are given in section 3. 

V. Conclusions 

POD and DMD methods can be utilised for analysing the LES data. The power of PODs fast convergence 
allows for the large scale structures to be isolated from the small scale structures in the turbulence. The 
DMD method has clear advantages over POD, as it strives for a representation of the dominant flow features 
with in a temporal orthogonal framework, while POD is based on a spatial orthogonal framework. This paper 
summerizes that POD and DMD methods will provide the experimentalist with solid tools, in quantifying 
important mechanisms in time resolved measurements of fluid dynamics. It is hoped that DMD and POD 
methods help in further understanding of fundamental fluid processes. In the future, effect of the time 
window between two frames of time resolved data on DMD method would be studied. Furthermore, the 
sensitivity analysis of the POD modes would be studied to show that the POD modes are independent of 
the number of snapshots. 
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